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Abstract 



The dielectric formalism is used to set up an approximate description of 
a spatially homogeneous weakly interacting Bose gas in the collision-less 
regime, which is both conserving and gap-less, and has coinciding poles of 
the single-particle Green's function and the density autocorrelation function 
in the Bose-condensed regime. The approximation takes into account the 
direct and the exchange interaction in a consistent way. The fulfillment of 
the generalized Ward identities related to the conservation of particle-number 
and the breaking of the gauge-symmetry is demonstrated. The dynamics at 
long wavelengths is considered in detail below and above the phase-transition, 
numerically and in certain limits also in analytical approximations. The ex- 
plicit form of the density autocorrelation function and the Green's function 
is exhibited and discussed. 
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I. INTRODUCTION 



The observation of Bose-Einstein condensation pj in trapped systems has given rise 
to an extensive experimental |2]|| and theoretical ||] study of weakly interacting dilute 
Bose-condensed gases. Recent measurements of the elementary excitations ]5]|| permit the 
detailed comparison with different theoretical approaches. Due to the existence of a num- 
ber of essential experimental parameters in the trapped systems, e.g. the trap frequencies, 
scattering length, size of the condensate and temperature, the validity of different approxi- 
mations depending on these conditions can be examined. Most of the theoretical approaches 
go back to the 1960s, and were originally derived for spatially homogeneous systems in view 
of their application to superfluid helium. For their application to the modern Bose-Einstein 
condensates in traps many of these approaches are now being extended to inhomogeneous 
systems and new approaches are also being developed. In the present paper we wish to 
explore such a new approximation scheme. It is of interest for the modern Bose-Einstein 
condensates, because measurements of quasi-homogeneous properties can also be made there 

0- 

From the theoretical point of view taking the homogeneous limit implies a great simplifi- 
cation. Due to the translation symmetry complicated integral equations reduce to algebraic 
equations; essential requirements like a gap-less single-particle spectrum (which applies only 
in the homogeneous system) and the particle number conservation can be expressed in the 
form of algebraic identities (Hugenholtz Pines theorem generalized Ward identities 
compressibility sum-rule etc.). In this form they provide important tests for the consis- 
tency of approximations obtained by selecting certain subsets of Feynman diagrams. Having 
defined an approximation or 'model' in this way analytical or numerical solution is possi- 
ble in the spatially homogeneous case without the need for further approximation and all 
branches of the collective modes - damped, strongly damped and even over-damped (purely 
damped) - can be exhibited. This is a great advantage over the corresponding calculations 
for the spatially inhomogeneous system, where further approximations are unavoidable and 
a similarly complete picture usually cannot be attained. 

The purpose of this paper is to examine the homogeneous limit of an approximation pre- 



viously derived for and applied to the inhomogeneous case [11 1. In contrast to the previous 



perturbative solution in the case of a harmonic trap [JI] we shall now avoid all further ap- 
proximations, after setting up the model by selecting the relevant class of Feynman diagrams, 
and derive and study all the excitation branches appearing in the homogeneous model. Fur- 
thermore, a number of exact theorems for the homogeneous system will be checked for their 
validity within our model approach. 

For example, in general the poles of the density autocorrelation function and the single- 
particle Green's function must coincide in the Bose-condensed region (see e.g. @). Phys- 
ically this is a consequence of the existence of the order-parameter describing the spon- 
taneously broken gauge-symmetry. In our approximation the required coincidence of the 
poles is a direct consequence of the derivation of the model within the dielectric formalism, 
which is designed to ensure just that. In contrast to a previous model studied within this 
formalism ( [13| , [12| and references therein) and therefore sharing this crucial property, the 



approximation or model studied here (and in |TT[ for the trapped system) takes into account 
the contributions of both the Hartree terms and the Fock (exchange) terms in the selected 
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graphs. Both types of terms are in fact of equal magnitude for the interaction via s-wave 
scattering in the experimentally realized dilute weakly interacting Bose condensates. 

In the next section we present the basic diagrammatic building blocks within the dielec- 
tric formalism as an approximation of the full system. The basic relations of the dielectric 
formalism are then used to obtain the final density autocorrelation and Green's functions 
determining the excitation spectrum. We also state there the Generalized Ward identities, 
which are exact relations ensuring the particle number conservation law in cases of sponta- 
neously broken symmetries. We prove in the appendix that these identities continue to hold 
in our model, which is a corner-stone of our model-building, because it provides an impor- 
tant consistency check. A further consistency check is the fulfillment of the compressibility 
sum-rule, which is demonstrated at the end of section II. In Section |T| we present the numer- 
ical results for the excitation branches at long wavelengths. Besides the modes we studied 
previously for the inhomogeneous system we find additional excitation branches belonging 
to excitations of mainly the thermal density. The behavior of these thermal branches will 
be studied in detail. For weak interaction they can be accurately approximated by simple 
analytic expressions, which is done in Section [TV]. For stronger interaction only numerical 
results are obtainable. They may suggest explanations for the appearance of the two dif- 
ferent branches measured at JILA JTJJ. Section [V] contains some conclusions and further 
remarks. 



II. THE MODEL 

A. The dielectric formalism 

Here we shall give some necessary background on the so-called 'dielectric formalism' at 
finite temperature, which was introduced for analyzing interacting Bose-condensed systems 
T5| , p!6|JT3| , p!7Hl9| , p^1 . An excellent exposition of the formalism and a complete list of references 



is given by Griffin \T2\. The presence of a Bose condensate breaks the gauge-symmetry of 
the system. As a consequence the density fluctuation spectra appear in the one particle 
excitation spectra and vice versa. In the dielectric formalism the corresponding coupling 
mechanism is explicitly exhibited. 

For sake of simplification we adopt the natural unit fi = 1 throughout the paper. In the 
Bose-condensed phase the usual decomposition of the Bose field operator ^ in a condensate 
wave function $ = an d the fluctuation operator $ with vanishing expectation value 

\j> = $ + <| (2.1) 
is made. (. . .) denotes the thermal averaging 

~ Try | e -/3(H- M 7V) 

(A) = : ^. 2.2 

The condensate wave function can in general be complex, but in the absence of vortices it 
can be chosen as real, without restriction of generality, and we shall consider this case in 
the following, for simplicity. The value of the chemical potential /i is fixed by the average 
number of particles in the system and is connected with the density of particles in the 



3 



condensate |$o| 2 by requiring ($) = 0. Imposing $ ^ we have the additional appearance 
of anomalous Green's functions and it is useful to introduce the matrix Green's function 

G a p(r,T;r',r') = -(T T [$ a (r, r)$J(r', r')] ) (2.3) 

with field operators in Matsubara representation $i(r, r) = $(r, r) and $2(7*, r ) = &(r, r). 
In addition to these one particle correlation functions we need the density autocorrelation 
function: 

XnnCr.TjrV) = -(r.fn^r^^r')]}. (2.4) 

Here n = n — (n) denotes the density deviation operator. The exact Green's function G a/ g 
can be obtained from Dyson's equation, i.e. by summing up a geometric series in terms of 
the Green's function G° a p of the noninteracting system and the self-energy E a( g: 

G a p = G° a p + G ary 'E ry gGgp. (2-5) 

The self-energy insertions are by definition irreducible i.e. cannot be split into two parts 
by cutting a single particle line. Different model approximations can be classified by the 
different kinds of interaction processes included in S a/3 . 

Beyond the definition of irreducible quantities it is equally useful to consider the proper 
contributions. Diagrams are called proper if they cannot be split into two parts, each con- 
nected to an external vertex or line, by cutting a single interaction line. They will be denoted 
by an additional tilde. The total density autocorrelation function x nn can be easily recovered 
from the proper contributions Xnn 

by just summing up a geometric series: 

Xnn(q,U) = - — r. 2.6 

1 - gXxm{q,u) 



Here u on the one hand denotes the Matsubara-frequency []2D| u = ioo n = z(2wr//3). On the 
other hand it is understood that the analytic continuation in u from the positive imaginary 
axis yields the retarded correlation functions, and therefore, if this is our purpose, uo will in 
the following always be assumed to lie in the upper half of the complex plane, g denotes 
the 2-particle T-matrix. In the whole temperature-domain of relevance here it can be taken 
as independent of q, to and be expressed by the s-wave scattering length a as 

Ana 

9 = 2.7 

m 

The importance of the use of the T-matrix in place of the bare interaction was first stressed 
in this context in Beliaev's work [21]]. (See |17]] for a discussion of ladder diagrams within 



the dielectric formalism at finite temperature.) As is well known since then the use of the 2- 
particle T-matrix in place of the bare interaction implies that a sum is implicitly performed 
over the ladder- diagrams of the two-particle interaction, and care has to be taken not to 
sum up the same class of diagrams again later-on, see e.g. [Q. 

In the dielectric formalism the approximations for E Q/ g and for \rm are related to each 
other via irreducible vertex functions A Q describing single processes of excitations out of the 
condensate and of relaxations into the condensate where the energy and momentum transfer 
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is described by an additional out-going interaction line. Unlike in the usual Green's function 
approach, where approximations are defined by the diagrams kept in the self-energies and 
Xnn as basic building blocks, in the dielectric formalism an approximation (or 'model') is 
defined by introducing building blocks for the proper and irreducible parts of the three 
quantities S a/ g, Xnn and A Q . In the following we call irreducible and proper quantities to be 
regular and denote them with an r in brackets. 

(r) (r) (r) 

Thus we start with specifying the building blocks Xn n , and A a and expressing all 



necessary quantities in terms of these three basic quantities, see e.g. ||12|| . We begin with 
the decomposition of the proper quantities G a p and x nn . Since a Green's function is proper 
if and only if all the self-energy insertions are proper we get 

G$(q,u) = (G%)-\q,u) - E«(, )W ) (2.8) 

with (G° /3 )~ 1 (q, fj) = 5 a/ 3 aua — (|^ — jj) where m is the mass of the atoms and we use 

the abbreviation a = 1 or -1 if the index a takes the values 1 or 2, respectively. Next we 
consider the contributions to x nn , which can be written in the form: 

Xnn = X^+^ ] G^4 } . (2.9) 

Here we have both, irreducible contributions x*m and contributions which contain at least 
one proper Green's function G a p. 

The expressions for x nn and S a/3 are obtained by summing up the corresponding geo- 
metric series: The result for x nn was already given in ( [2.6|) , and we can immediately turn 
to the elements of £ a( g, which we decompose into their proper and their improper parts as 
follows: 

Z aP (q,u) = ^(q,u)+A£\q,u) -A<g\q,u;). (2.10) 

($) = implies that £j (0, 0) = 0, where (q, oj) is the sum of the regular diagrams 
with one external line (tadpole-diagrams). 

In addition the conservation laws and corresponding sum rules must be woven into the 
formalism by ensuring that the approximate expressions for the regular parts still satisfy 
the exact generalized Ward identities. The identities we want to consider in the following 
are a consequence of particle number conservation which becomes nontrivial because of the 
broken gauge-symmetry. Within the dielectric formalism the corresponding Ward-identities 
were first introduced by Wong and Gould || at temperature T = and later developed 



further for finite temperature (see e.g. ||T8|,|l2"f and further references given there) and used 



as exact relations between the building-blocks of the dielectric formalism to be satisfied by 
any consistent approximation. They can be obtained by inserting the continuity equation 

^ = -^j(,,r) (2.11) 
or m 

into the expressions for various time-ordered correlation functions 0, [12] and read: 
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ujt&\ 



q,u) 



m 



uj — a( /i) 

2m 



- \/n. 



UJ] 



"X&W) = j^X&W) - V^c [AS r) ( 9)W ) - 4\q,u] 
uXjni<li u) = —[Xjj\q, <*>) + m(n c + n th )] - v 7 ^ f A i (r) (9) w ) - ^(Q, 



uj 



m 



(2.12) 

(2.13) 
(2.14) 



The demonstration that eqs( [2. 12 - 2.14 ) are satisfied by our model approximation is crucial 
for ensuring its consistency and given in the appendix. 

A further important requirement for the consistency of any approximation is the com- 
pressibility sum-rule [jTOf 



lim Xnn(&,0) 
|fc|— >0 



dn t 



[2.15) 



To these exact requirements should in our opinion be added a further one, which is a con- 
sequence of Galilei-invariance and will be discussed in the final section. 



B. The Model 

We shall now introduce and motivate the particular approximation (or 'model') within 
the framework of the dielectric formalism we wish to study in this paper. Our aim here is to 
give a self-consistent Hartree-Fock theory of the Bose-gas, not a perturbation theory in the 
coupling constant to any given order. This aim seems legitimate in view of the success of the 
Hartree-Fock approximation in many-body physics in general, its limitations as a mean-field 
theory notwithstanding. 

Let us first consider the Bose-gas in the non-condensed phase. Previously this system 
was studied in the framework of the dielectric formalism within the Hartree-approximation 
HI , which corresponds to the self-consistent theory summing all the bubble- diagrams, but 



neglecting exchange-contributions in the thermal cloud.. Here we would like to go a step 
further and study the self-consistent Hartree-Fock theory summing the bubble-diagrams and 
the exchange-bubble diagrams. For the uncondensed phase this approximation is, of course, 
standard, both for the single-particle Green's function, where it corresponds to the standard 
Hartree-Fock approximation, and for the density response-function, where the summing of 
the exchange-bubbles is discussed e.g. in [|(J. 

The dynamics of the thermal cloud is described by the Hartree-Fock Green's function, 
which takes the form 



with 



2m 



With these Green's functions the lowest order contributions to the regular density auto- 
correlation function Xnn (see fig.[3]) are given by the so-called bubble-diagram 
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x l(q,uj) = -^j±y 3 J tPpGFtokjJGg^q-au-iuJ. (2.18) 

For the sake of better readability we will use the abbreviations p = (p, iu n ) representing the 
momentum p and the Matsubara frequency iu n . 

In order to treat the exchange-contribution of the two-particle scattering in Xnn on the 
same footing as the direct contribution we have to sum up the diagrams for Xnn as in figfj], 
which define a simple geometric series. Since the interaction g is given here just by a constant 
the result can be written in terms of the expression ( 2.18 ) in the form 



Y° (q) 

Xnn(q) = = o i v ( 2 - 19 ) 



When this result for Xnn is substituted in eq. (|2.6| ) it yields simply 

i.e. a doubling of the effective coupling constant in the denominator. Concerning this 
treatment one may legitimately ask, why the Hartree-Fock Green's functions rather than 
the free-particle Green's functions are used for the internal propagators in the bubble and 
exchange-bubble diagrams. The reason is one of consistency. Without the use of Hartree- 
Fock propagators one could not achieve consistency between statics and dynamics, which 
is incorporated in the compressibility sum-rule (|2.15|) (nor, as we may add, the consistency 
with Galilei-invariance, as we discuss in the final section). 

Let us now turn to the description of our approximation, or 'model', also for the con- 
densed phase. It is defined as the natural extension of the Hartree-Fock model we use in 
the uncondensed phase. Thus, for the internal lines we use again the Hartree-Fock propa- 
gators. In this way we again achieve consistency with the compressibility sum-rule as will 
be shown later. Furthermore, we again sum up the bubble and exchange-bubble diagrams 
for all proper quantities, which now also include the proper vertex function A^. Therefore 
we find 

y( r )( a ) = *nn(g) (2 21) 



which is almost the same expression as found in the uncondensed phase, see also |27[| . 

For the proper vertex function A„ ■* (q) the zero-loop diagram is given by the trivial vertex 
function A° = $o- By the same reasoning as for \nn we obtain a further geometric series: 

AM(8) = T^W" (2 - 22) 



illustrated in Fig.^. 

The regular self-energy contributions given by Ti^l(q) = E°'g (g) + ^^' ) (g) are presented 



in Fig.|], where 

C ) = (^l $ o| 2 + 2^ th )(j°) • (2.23) 

\ ' a.B 
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a/3 



We now simply insert the expressions (ggg) , ( gjg) , (gjg ) and (|2T21| ) in Eq. ( ^T0|) 



and obtain straightforwardly the expressions for E Q/ g(g) of our approximation, which are 
not written out here to conserve space. Using our results for En and E 12 it can be easily 
checked that the Hugenholtz-Pines theorem fi = E n (0) — E 12 (0) is satisfied. 

The tadpole diagrams are given in fig.£|. The condition E^ = leads to a relation which 
can be solved for the chemical potential 

/i = gn c + 2gn th + /x (0) , (2.25) 

where we added the chemical potential of the ideal Bose-gas, which vanishes in the 
Bose-condensed phase, in order to make the expression valid also in the uncondensed phase, 
and where we define 

"* = "£ / j^G HF (j>,iu n ). (2.26) 

That n t h is consistently interpreted as the thermal density of the particles in our approxi- 
mation will be shown by checking the Ward-identities and the compressibility sum-rule. 

One might raise the question why in the diagrams for E^ and the diagonal elements Eq2 
exchange-ladder terms with two condensate lines at the same interaction line are not also 
included. The reason is that because of our use of Hartree-Fock propagators in the internal 
lines these exchange-diagrams are already contained in the first-order proper exchange- 
diagram. As a result the chemical potential fi in eq.( [2.25| ) is simply given by the lowest 



order terms, but of course these must be evaluated with the Hartree-Fock propagators. 

Using the building-blocks we have specified it is now straight-forward to obtain the 
explicit expressions for the Green's function, namely 



G n (k,u) 



(w + £ - a*<°>)(1 - 2g X ° nn (k, w)) + gn c (l + 2g^ m (k, w)) 



A(fc, u) 

with the denominator 

A(fc, u) = (u 2 - (^ - ^) 2 ) (1 - 2g X Uk, «)) - ^n c ^(l + 2g X ° nn (k, w )). (2.28) 
The density autocorrelation function can be given similarly as 

, (w 2 -(f -ii(°)) 2 )xi (k,uj)+2n c ^(l + gx° nn (k,uj)) , x 

Xnn (k, U) = Dn \ , r ^ . 2.29 

A(k, to) 

The poles of the Green's function and the density autocorrelation function are given by 



A(k,uj) = 0. In the uncondensed phase we can put n c = and the result (|2.29 ) reduces 
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to eg. ( |2.20| ) with the poles given by gx° n (fc,co>) = 1/2. These dispersion-relations will be 
analyzed in section |T|. 

It is manifest from our explicit expressions that the poles of the Green's function in the 
Bose-condensed phase are indeed the same as those of the density autocorrelation function. 
On the other hand, because of n c = the poles of Xnn(k,u) and Gn(k,u) become different 
in the uncondensed phase: the factor (u + |j- — — 2gx° n (fc, (J)) then cancels from 

the numerator and denominator of Gn(fe,w) and the remaining single-particle poles at 
to = [k 2 /2m) — are those of the free Bose-gas. 

It is important to remark here, that summing up the bubble and exchange-bubble dia- 
grams, after a suitable rearrangement of the resulting expressions, has not lead us to terms 
of higher than the first loop order in the numerator and denominator, both in the condensed 
and the uncondensed phase, if note is made of the fact that the quantity gn c is of zeroth 
loop-order in the condensed phase. However, our treatment is not a perturbation theory 
in the coupling g, i.e. loop-order is not the same as order in the coupling g, because a 
self-consistently determined propagator, namely the Hartree-Fock propagator, is used in the 
bubble-diagram defining x° n , which depends itself on g. Rather, we claim that the theory 
given here is the self-consistent Hartree-Fock theory also in the condensed phase. This claim 
is based on the fact that the theory satisfies all consistency checks, namely the identity of the 
poles of the Green's- and density response-function, the Ward-identities, which also contain 
the f-sum rule, the compressibility sum-rule, and of course the Hugenholtz-Pines theorem. 
In addition it satisfies a consistency check derived from Galilei-invariance, as discussed in 
the final section. It should be noted that, apart from the Hugenholtz-Pines theorem, none 
of these consistency checks are satisfied in the usual Popov approximation. (But it is inter- 
esting to remark that they can also be shown to be satisfied in the simpler Hartree-model 
TB|JKj|| within the dielectric formalism, for which our treatment therefore provides the logical 



consistent extension which includes exchange). 

Surprisingly for us the comparison of our result ( J2.29 ) for the density autocorrelation 



function with corresponding results given by Minguzzi and Tosi |28[ yields complete agree- 
ment. This is nontrivial, since the theories are formulated in quite different ways and arrive 
at the result for the density autocorrelation function in very different manners and after 
considerable rearrangements. It is interesting to learn from this agreement that the theory 
formulated by Minguzzi and Tosi [28| has a foundation within the dielectric formalism and 



permits a description in terms of the Feynman diagrams from which we derived our results. 
This is a strong hint that there is just one consistent Hartree-Fock theory not only in the 
uncondensed but also in the condensed phase. 

An obvious advantage of the dielectric formalism over the linear response theory used 



in [28| is the fact that besides the density response-function the Green's function is also 
obtained explicitly as displayed above, which would be outside the scope of the theory given 
in PBfl . With the Green's function our theory gives explicit expressions for autocorrelation 
functions of amplitude and phase of the order-parameter. The former equals — G + , the latter 
— G_, with 



G + (k,u)= Gn(k,u) + G 12 {k,LU 



^ + £-^ (0) )(l-2 ^n(fc^)) (23Q) 
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G-(k,u)=G 11 (k.u)-G vl (k.u) 



(uj + ^- /iW) (1 - 2gxUfc, u)) + 2gn c (l + 2^° n (fc, gQ) 

We still need to evaluate the thermal density n th , the chemical potential \i and n c . 
By integrating the Bose factor 

over the momenta we get n th = g 3 / 2 (z)/\ : l h , where g3/ 2 (z) is the Bose function [PD| , 2; is the 
fugacity, defined in terms of an effective chemical potential p = p — 2gn tot via z = exp(Pp), 
and Ath = \/27r /rahsT denotes the thermal wavelength. In the uncondensed region the 
effective chemical potential p coincides with the chemical potential /i*- ** of the ideal Bose- 
gas. For the condensed phase, using \i = g n c + 2gn t h, the effective chemical potential is 
p = \i — 2g (nth + n c) = —gn c and z = e~^ nc . For a given total density n to t = n c + n t h the 
equation of state in the condensed phase is therefore given by the implicit equation, which 
is equivalent to the equation of state derived first by Huang et al. 



n to t=n c + ^g 3/2 (e-^). (2.33) 
A th 

It can be seen from this equation that within our approximation for G HF , the phase 
transition is no longer continuous^. : E.g. for T = T c , where T c is defined as the critical 
temperature of the free Bose-gas for the same particle-density by n tot = 5f 3 / 2 (l)/Aj s h (T c ), 
there are two solutions n c = and n c = An c , with 

47] _g(\ c h )~ 3 ,, c s 3 

An c (T c ) (A t c h ) 3 = ;^ AC) -3 N9 + 0{{^f^fl\ (2.34) 

( l_ C( l /2 )^_)2 k B T c 

A^ h denotes the thermal wavelength at the critical temperature. Even for comparably weak 
interaction gn tot /kBT c = 0.1 this jump already amounts to 16% of the total density, so it 
cannot be neglected by any means. The coexisting values of the total density for given 
values of the chemical potential in the two-phase region have to be determined from the 
Maxwell-construction and will be discussed below. 

Let us first consider the correlation functions in the static limit, i.e. at u = 0, in more 
detail. We find with or without condensate 



G.&OH-^J— . (2.35) 

2m r 1 



With a condensate //°) = and this correlation function displays the usual infrared singu- 
larity ~ k~ 2 associated to the spontaneously broken gauge-symmetry. For the amplitude 



i A discontinuity is also found in the familiar Popov-approximation [30,|31|, see e.g. [32|,33,22,34| 
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autocorrelation function at w = we obtain in the uncondensed and the condensed phase, 
respectively, and displaying also the asymptotics for \k\ — > 0, 



2rn 



(2.36) 



G(k0 )= l-2gX°n(fc,0) , -1 T <T 

r+ fc2 (l-2^ n (fc,0)) + 2^ c (l + 2^ n (fc,0)) 2Eia(fc,0) * ' 



2m' 



with 



l-2^ n (fc,0) 



S 12 (fc,0) gn c (l + 2g X nn (k,0)) 



[2.37) 



We examine the stability of the model below T c . As can be seen from the results for 
G + (k, 0) an instability may occur if for small \k\ — ► we satisfy gx^nik, 0) = —1/2. In order 
to see where this condition for instability is satisfied we evaluate Xnn(O'O) = ^ m Xnn(^'O) 

from Eq. ( |2.18|) , finding 



^°„(0,0) = -^ 1/2 (e-^). (2.38) 
A th 



This is a negative quantity, which, for sufficiently low temperature, lies in the interval 
[—1/2,0], but which moves monotonously in the direction of the instability at -1/2 with 
decreasing condensate density for fixed temperature. The instability condition 5X^(0, 0) = 
— 1/2 is equivalent to the condition that for fixed temperature the chemical potential ii 
exhibits a minimum ^( min ) as a function of n c . In Fig. |5] we plot this functional behavior. 
The fact that the instability condition, which we have obtained here from the dynamics 
described by G a p, fits completely with the form of the equation of state ( ^.33[ ) underscores 
the consistency between statics and dynamics in our model. 

In Fig. H we plot the branch corresponding to the chemical potential of the uncondensed 
Bose gas in Hartree-Fock approximation together with the branch for the condensed system. 
We only get solutions for the uncondensed Bose gas if fx = //^ < or equivalently \i < 
^(max) _ 2gA^ 3 C(3/2). Combining these results due to the first order transition we find a 
region of mult i- stability between the condensed and the uncondensed Bose gas for chemical 
potential values /i( min ) < fi < /i^ max \ In this region there are actually two possible non- 
vanishing values for the condensate density n c ^ for given temperature and chemical 
potential. But due to the stability condition gx° n (0,0) > —1/2 only one of these solutions 
proves to be stable. The coexisting values of n tot at a given temperature are fixed by 
the Maxwell-construction for the /^(n tot )-curve. The resulting value fi coex . of the chemical 
potential then also determines the value of n c coexisting with n c = from Fig. [5|. 

The equation of state as obtained from Eq. Q2.33Q and plotted in Fig.|| is fully consistent 
with the density- density correlation function ( |2.29| ): the compressibility sum-rule ( |2.15| ) 
makes this connection between statics and dynamics. With our present result for x„„(k, 0) 
it takes the form 

dn tot _ 1 1 + gxLM (23Q) 
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That the derivative {dntot/dfi)T can indeed be written in this form can be checked by a 
short calculation using the identity ^#3/2 (e^ 3 ) = — /3<7i/ 2 (e^) and eqs. (|2.33|J2.38| ). Thus, 
from the right-hand side of equation ( |2.39|) , which is obtained from the dynamics, a static 
property on its left-hand side, namely the equation of state, can be re-derived by integration 
with respect to ll. 

The slope ^-^\t for n to t — > 00 equals the coupling constant g according to both sides of 
Eq. ( gggD - The border of instability ^\ T = occurs where the density autocorrelation 
function has its pole, i.e. for 2gx° n (0,0) = —1. It should be noted that this is the same 
point where the longitudinal order-parameter response-function G + (k, 0) has its singularity, 
and also where (dn c /d/j,)T becomes singular, as follows from the relation dfi/dn c = g(l + 
2<7x° n (0, 0)), which can be easily proven. The infinite slope ^-^\t = 00 of the unstable 
branch occurs where Xnn(0,0) vanishes, i.e. for gx^ n (0,0) = —1. 



III. DYNAMICS AT LONG WAVELENGTHS 

Now we want to investigate the dynamics at long wavelengths. We consider the limit 
to — > and \k\ ^ and introduce the complex velocity of sound by oj/\k\ = c. The 
dispersion-relation in the condensed phase, given by the poles of the density autocorrelation 
function and the Green's function, takes the form 

(c 2 -c 2 B ) = (c 2 + c 2 B )2g X ° nQ (k,u;), (3.1) 



where we introduced the Bogoliubov speed of sound cb = ygn c /m. Eq. ( |3.1| ) will be 
the central equation we shall analyze in this section. The response function (|2.18 ) for a 
homogeneous system can be rewritten in the well-known form 

(h , f d 3 p fo(E HF (p))-f (E HF (p+k)) 
XmA ' ' J (2tt) 3 u-(E HF (p + k)-E HF (p)) ' 1 ' 

where fo(E HF ) = (exp(/3(E HF )) — 1) _1 is the Bose-distribution. For T < T c we have 
E HF {k) = k 2 /2m + gn c and p, = —gn c ,. 

The dispersion-relation for the density-fluctuations in the uncondensed phase is formally 
also contained in Eq. ( |3.1| ) putting there cb = and has one branch with c = 0, corresponding 
to the fact that the single-particle dispersion-law in the uncondensed phase is proportional 
to k 2 , and another branch given by 

l = 2g X ° nn (k,uj). (3.3) 

The response function ( |3.2j ) in that region is evaluated with the effective fugacity z = 
, according to p, = [l — 2gn tot = , where the chemical potential jx is determined 

thermodynamically as described above. 

For long wavelengths, i.e. for wave vectors with |fc|A t h < 1, we can approximate in the 

denominator of the response function fl3j) E HF (p + k) - E HF (p) w k-p/m. The mean 



field interactions for our two models cancel in this difference. The difference of the Bose 
factors in the numerator can be be approximated by a gradient. To evaluate ( |3.2|) after 
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these approximations further we choose k in x-direction and integrate first in p y ,p z with the 
result 



= ^ j *. W | fc | p : Pi/m M^ f fe))- 

The response function only depends on the ratio oo/\k\ , which defines the (generally complex) 
'speed of sound' defined by the ratio c = uj/\k\ in the long wavelength limit. So equation 
TO) is an implicit equation for c . Scaling p x by \/2mkBT and measuring c in units of the 



thermal velocity ct = ^/2ksT /m , the response function reads 



Xl(c) = xL(c/c T ) (3.4) 
k B T 



with the dimensionless response function 



X° nn (s) = -L I dt— , (3.5) 

- 'n s - t e 



defined in the upper half of the complex s-plane, where the variable s denotes the quotient 
s = a;/(cy|fe|). The temperature dependence resides in the prefactor in ( |3.4j ) and in the 
fugacity z. 

The integral can be evaluated by the methods of residua, see [pj , and we refer to this 



paper for further details of the calculation. The result is the expression for the dimensionless 
response function 



SZ 1 7T 

Xnn(s) = -9l/2(z) ~ \f*i — + J- 

e s — z y 7 



7 s 2 + 7 



with poles at 



— - — + S — . (3.6) 
*-^n=i s z - a z n a n s z - b z n b n 



b 



ii 



and with 7 = = \/3fi\ = — ln(z) , and <p n = arctan(27m/7) for n > 0. In the Hartree- 
model |T^] 7 was in the condensed phase, while here it is non-zero in the condensed and 



uncondensed phase. This may seem like a small difference, but his would be misleading, 
because a non-zero value of 7, even if it is small, leads to qualitative and by no means 
small differences for the response-function in the limit k — >• 0. For small 7 one can expand 
<?i/ 2 (e^) = vWt+C(1/2)+0( T ). 

When we ask for solutions of (|3.1| , |3.3|) , we are looking for eigen-modes of the system, 
which, for physical reasons, have to decay (rather than grow) exponentially and the poles 
we look for have therefore to be located in the lower complex half-plane. In Eq. ( |3.1| ) the 
analytical continuation of the integral (|3.5| ) from the upper to the lower complex half-plane 
of s has therefore to be used. 

With this explicit representation of the response function we are now able to determine 



the solutions of ( [3.1|) and of (|3.3| ) below and above the phase transition numerically. Since 
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the terms in the sum in (|3.6| ) decay only as n~ 3 ^ 2 , the sum converges very slowly and it 
is not possible to cut the sum at some large finite n cut . Instead we can use a continuum 
approximation integrating over all terms in the range n G [n cnt , oo[ . As a numerical check we 
have reproduced all our data with the alternative representation of the response function by 
performing integral ( |3.5|) numerically with s in the lower complex half-plane and subtracting 
the term for the analytical continuation: 

xL( s ) = ZT72 / dt ~~ + / ~ J Z .. • ( 3 ' 7 ) 



n 1 '* s — te"~—z e s 



The numerical solutions are plotted in Figs.|10L 11}. These figures show the real- and imaginary 



parts of c for various branches of damped modes as a function of temperature above and 
below the phase-transition, and how these branches bifurcate as the poles corresponding to 
these modes move in the complex plane. For the weakest interaction strength, shown in 



Fig.|10|, various approximations are possible to permit an analytical understanding of most 
of the structure shown. This will be described in the following section. Here we discuss the 
numerical results shown in the figures. We use the condensation temperature k^T c of the 
ideal Bose-g clS cLS cL convenient energy-scale near which the phase-transition occurs. 

Let us begin with the case of weak interaction with gn tot /kBT c = 0.1 for our two models, 



shown in Fig.|10|. At high temperatures T > T c we typically get complex solutions for c with 



|c| ^> ct , and with finite real and imaginary part. Then, at a certain temperature T , which 



is above T c for very weak interaction but moves below T c for stronger ones as seen in Fig. [1 I 
the real part of the velocity of sound vanishes, the mode becomes over-damped below T c , and 
the imaginary part bifurcates into two different branches, as seen in Fig. |10|. Density-waves 
with a certain wave vector no longer propagate, but decay as a mere relaxation with two 
different decay-rates describing a short-time decay and a long-time decay. Both damping- 
rates are proportional to the wave vector. An analytical understanding of the bifurcation 
based on a suitable approximation to the response function for large arguments \s\ 3> 1 will 
be provided in the next section. 



The finite jump ( |2.34j ) of the condensate density at the first order transition in the model 



results in discontinuities of the two purely imaginary thermal branches. The discontinuity 
in the lower branch is visible in Fig. |It], while in the upper branch it is very small and not 
discernible. Even at lower temperatures this upper branch hardly deviates from the result 
above T c as can be seen in Fig. |TU[ The reason for this is that for weak interaction ( fl.2| ) also 
implies |c| 3> cb , and the dispersion-relation ( |3.1| ) below T c reduces to the equation ( p.3| ) 
above T c . So c from (fO|) turns out to be also a solution of (|3.1|) below T c and can be seen 



in Fig. [10] as the solution independent of the temperature for the whole condensed region. 

Below the phase-transition a new propagating branch with non-vanishing real part of c 
appears, which is the Bogoliubov mode and, as discussed in the next section, has a velocity 



close to Cb, as is also visible in Fig. 10. In the units chosen in the figures this Bogoliubov 



branch for T — > converges to the value \J gn tot /k B T c . 

Let us now turn to the case of slightly stronger interaction in Fig. 11]. The most 
remarkable new feature compared to the very weakly interacting case is, that the bifurcation 
of the thermal branch occurs now at some temperature T below the phase transition. We 
can see the real parts of the propagating thermal branch and the Bogoliubov branch crossing 
at some temperature between T and T c . Though the mode described by the thermal branch 
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is still propagating it is strongly damped, its real part being smaller than its imaginary part. 
The bifurcation of the Bogoliubov branch near T c does not occur anymore, the Bogoliubov 
sound remains propagating at the phase transition. Again the smaller of the two damping 
rates which have bifurcated from the thermal branch is still larger than the damping of 
the Bogoliubov sound. However, the discontinuity in the condensate density at the phase 
transition already amounts to 40% of the total density, so this behavior near T c has only a 
restricted physical meaning. 



IV. ANALYTICAL SOLUTION OF THE DISPERSION RELATION FOR WEAK 

INTERACTION 

A. Thermal branch 



The bifurcation shown in Fig.|10] can be understood by the following approximation to 
the response function: for large arguments \s\ ^> 1 we can expand the integrand in ( |3.7| ) in 
powers of (t/s) 2 and get as the dominant behavior for large \s\ 

,o„ w ^_ 2 ^_££_ + (i). (4,) 

For weak interaction g n tot <C k B T and for T > T c the real part of s is much smaller than the 
imaginary part, as seen numerically, and the full response function can be approximated by 
the second term in ( |4.1| ) only. This approximate response function, taken as a function of a 
purely imaginary argument, is a real convex function, and only if condition ( ^.3| ) is fulfilled at 
its minimum, a purely imaginary solution is possible. From the temperature dependence of 
the chemical potential above T c it follows, that this is possible only below some temperature 
To, which one can determine as Tq/T c = 3.37 for gn^tlksTc = 0.1 , for example. Above To 
the approximate response function gives rise to the characteristic growth of the real part of 
c, see Fig. [10]. In the limit \s\ > 1 we even can approximate Xnn( s ) ~ 2y/nis, and the 
simple solution independent of the temperature follows as 

i k B T m 

c = -^^5 CT = ~ a - (4 ' 2) 

This can be seen as the upper branch of the bifurcated imaginary part for small temperatures 
in Fig. [TTJ. The lower branch decreases linearly with temperature and finally enters the 



opposite region with \s\ = \c/ct\ C 1.. In this region the behavior can be described by a 



different approximation of the response function [TBI]. The first three leading terms of 



in s, 7 ordered by magnitude can be summarized as 



,~,o ,„\ _ h X \ , , ^i„i2 



— , +v ? <( 5 ) + V +oa 'i')- (4 ' 3) 

From the dominant, first two terms the approximate solution in the limiting case |c| ct 
above the phase transition follows as 
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Except at high temperatures this expression agrees very well with the numerically deter- 
mined data in Figs.^. The behavior linear in T — T c is given by the first term in the numer- 
ator, since the chemical potential can be expanded as —\i = C(3/2) 2 9/cb(T — T c ) 2 /16nT c + 
0((T — T c ) 3 ), the offset is given by the second term in the numerator. 

Finally we note, that the bifurcations discussed above can also be followed analytically 
if use is made of the approximation (|4.3|) for the density response function. 



B. Bogoliubov-branch 

Let us turn to the region below T c and show that in both of our models we have a 
branch with c very close to eg. In the limit |c| <C ct we can use again the approximate 
response function ( |4.3| ). According to the Bogoliubov theory at vanishing temperature sound 
propagates with the Bogoliubov speed c B . Using this in (|4.3|) we see, that gx^nips) is of 
the order 0{sb ) 1 1 ^ 2 ) <C 1 with sb = cb/ct = ylj^. So with gXnn( c B) <C 1 it follows 
from ( |3.1|) , that one solution is near the Bogoliubov speed of sound c# . The corresponding 
numerically determined solutions can be seen in Fig.pH| 

Being close to cb the Bogoliubov-branch can be determined analytically from pertur- 
bation theory. To this purpose we consider the right hand side in ( |3.1| ) with the response 
function ( |4.3|) as a perturbation to the Bogoliubov speed of sound c = cb ■ We get the 
complex correction 

c-c B = 2c B ^(^^-C(i) + ^p). (4.5) 
k B T \s B + ty^f 2 2 / 

The first term gives the leading order in g and is proportional to the temperature. The 
damping rate of Bogoliubov sound waves with wave vector q due to this term is 

T 1 = ^k B Taq 1 (4.6) 

with the s-wave scattering length a from g = 4ira/m. This damping rate is of the same 
order as the damping previously determined in [13|], which only differs in the numerical 
prefactor of one instead of 4/3 , and also agrees approximately with the damping rate in the 
Beliaev approximation extended to finite temperatures by Shi and Griffin |22| , but obtained 
also by other methods in the intermediate temperature region p3|-p6[]. The difference is 



that in ( |4.6|) the prefactor 4/3 should be replaced by 3ir/8. For the frequency shift we 
obtain Au = — y/2T, a result which is in good agreement with a result of Fedichev and 
Shlyapnikov p5| , who obtain Au/T = — (28/37r 3//2 ) = 1.67..., and also of Giorgini |35|], who 



gets Aoj/T ~ —1.8 This shift, which is negative and linear in the temperature, reduces the 
speed of sound compared to the Bogoliubov approximation. These results are in agreement 



with the measurements of temperature dependent frequency shifts of discrete modes |]14 



which were found to have negative sign for the m = 2 mode, and also for the m = mode at 
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intermediate temperatures. Using the full response function and taking the inhomogeneity 
of the system into account an explanation of these frequency shifts and the damping rates 
was already given in 111 . 



V. DISCUSSION AND CONCLUSION 

Let us briefly summarize the results of this paper and then draw some further conclusions. 
We have presented here within the framework of the dielectric formalism a consistent micro- 
scopic model of the weakly interacting Bose-gas including exchange. We have shown that 
a consistent treatment of exchange processes is achieved by using Hartree-Fock propagators 
for the internal lines of diagrams and summing up the same classes of diagrams for different 
quantities. As far as the density correlation function is concerned we obtain results which, on 



the general level, are equivalent to earlier reults by Minguzzi and Tosi [28]. This agreement is 



nontrivial, because our starting point is quite different from their's. However our treatment 



is more general than that in p8| because it also gives the single-particle Green's functions, 
which by construction have the same poles as the density- correlation function, displaying a 
gapless single-particle spectrum, i.e. satisfying the Hugenholtz-Pines theorem. The agree- 



ment of our result for the density correlation with that in EH] uncovers the diagrammatical 



basis of the equations written down there and therefore opens up the possibility for future 
systematic improvements. The rational basis of our 'model' or approximation is its consis- 
tency with general requirements, which are very nontrivial to be satisfied simultaneously. 
Thus we demonstrated explicitely that the compressibility sum-rule is satisfied, ensuring 
the consistency between statics and dynamics of the model; and the Ward-identities, were 
checked, ensuring the consistency between particle number conservation (and the f sum-rule) 
and the spontaneously broken gauge-symmetry. 

To these consistency checks, which were in detail discussed, a further one may be added 
which we have not yet discussed, but which we deem to be of no less importance, because 
it derives from a further symmetry of the system - Galilei-invariance. Galilei-invariance is 
most easily considered in a spatially confined system, because it then simply implies the free 
motion of the center of mass, if the confined system as a whole moves. In Bose-condensed 
systems spatial confinement is naturally achieved by imposing an external, but spatially 
fixed trapping potential. The system can then not move as a whole, but can still move in 
the external potential. The motion of the center of mass is then no longer free, and in general 
it is not even separable from the other degrees of freedom in quantum mechanics. However 
in the special case of an external harmonic potential the center of mass motion is separable 
and is simply a harmonic oscillation in the external potential. This fact is the content of 



the Kohn-theorem ||36|| . In the limit where the spring-constants of the external harmonic 
potential are sent to zero the harmonic oscillation of the center of mass tends to the free 
motion required by Galilei-invariance. Therefore, if the system satisfies the Kohn-theorem 
in a fixed external harmonic potential (possibly with infinitesimal spring-constants), Galilei- 
invariance is ensured if the external potential is switched off. We wish to point out here, 
that in addititon to the other consistency-checks already discussed, also this check is passed 
by the approximate model discussed in the present paper. This has been shown in |[37f , 



where the fulfillment of the Kohn-theorem within the approximate model was explicitely 
demonstrated. It can be seen from the proof given there that this test is quite sensitive and 
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would e.g. fail if. another definition of nth would be chosen, e.g. by using in the definition the 
full Green's function rather than the Hartree-Fock propagators, or if propagators different 
from the Hartree-Fock propagators would be used in the internal lines. 

Let us briefly discuss also the limits of validity of our treatment. The theory we have 
given is a mean-field theory and therefore restricted to the collisionless domain ujt co u ^> 1, 
where t co u ^> 1 is the mean collisiontime. It is also restricted to sufficiently high temperature 
ksT > gn c , because in the low-temperature domain ksT < gn c scattering processes at wave- 
number smaller than the inverse Bogoliubov coherence length £g = (87rn c a) -1 / 2 make an 
important contribution, and our use of the Hartree-Fock propagator for the internal lines 
would lead to a qualitatively wrong temperature-dependence. The mean-field character of 
our theory also prohibits its application close to the phase transition, which occurs at a 
temperature close to T c for the weakly interacting Bose-gas. The Ginzburg-criterion for the 
validity of a mean-field description here takes the form [25| \T — T c \/T c > (n^a 3 ) 1 / 3 . Indeed, 



our model, if extrapolated to temperatures near T c , would predict a first-order transition 
within the transition-region specified by this criterion, but since this is clearly outside the 
limit of validity, it is of course not a prediction of the model. For some purposes, like 
following the fate of the various excitation branches as the phase-transition temperature 
is crossed, it would certainly be nice to have also a mean-field model including exchange 
and satisfying all the consistency checks we have discussed and also giving a second-order 
mean-field transition at a critical temperature near that of the ideal Bose-gas. This goal, 
however, is not met by the approximation we have discussed here, and further work may be 
required to achieve it eventually. 

In summary, the results obtained here identify the model we introduce as a rather sat- 
isfactory while still manageable microscopic description of a weakly interacting Bose-gas in 
the collisionless regime, except at very low temperatures and very near to T c . 

Besides checking in detail the consistency of our approximations, we have presented and 
discussed a detailed numerical and partially also analytical study of the dispersion relation of 
the joint single-particle and density fluctuation modes below and the separate single-particle 
and density fluctuation modes above T c . The results for the complex ratio c = u)/\k\ have 



been summarized in Figs. |10yiT| . The dispersion relation found depends qualitatively on the 
strength of the coupling. If the latter is weak like in Fig.pl], then purely damped modes exist 
from a region above T c down to the low-temperature regime, besides the propagating and 
weakly damped Bogoliubov-mode which exists only below T c . For stronger coupling, like in 
Fig.|TT|, there is a propagating, damped mode from above T c down to a finite temperature 
T somewhat below T c . Only below T this mode also becomes purely overdamped, like for 
the weak-coupling case. The Bogoliubov mode in the condensed phase exists also for strong 
coupling, only with higher frequency and larger damping. 

Though the main purpose of the present work has been a theory for the homogeneous 
Bose gas it is interesting to contrast the results presented here with the measurements of 
and theoretical results for the temperature dependence of discrete frequencies for trapped 
condensates [H], [3S|, The results for the real part of the velocity of sound in Fig.|TT| 
are similar to measurements of the frequency shifts and damping-rates for the m = and 
m = 2 modes at intermediate temperatures ]14|]. The measured frequency of the m = 
mode was found to first decrease with increasing temperatures for low temperatures, but 
then it suddenly increased again at higher temperatures T > 0.6 T crit . It was suggested, that 
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the increase might be due to the crossing of the Bogoliubov mode with another mode of the 
l38| , [[391 . Such a second mode was actually found in |38 |, J39| using a 



thermal cloud 14 



solution of the kinetic equations. Here we also found a second branch of the velocity of sound 
for the homogeneous system, see Fig.[y], but it is found to be strongly damped, in qualitative 
distinction from what is found in the trapped system. However, a measurement along the 
lines of , testing the local properties of the Bose-gas, could actually check our results for 
the thermal branch of the collective modes in the homogeneous system. In fact it was already 
remarked in |7j that above T c there was no clear evidence for a propagating sound-wave, 
which is in qualitative agreement with the non-propagating nature of the collective mode 
above T c which we have found in the present paper. 
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APPENDIX: 

We wish to show here that the identities ( |2.12| - PTT4"D are indeed fulfilled in our model 
approximations. The specific choice of the building blocks for the vertex functions A^ and 
Aa simplifies the structure of these identities considerably. First the second term on the 

(r) (r) 

right hand side of the second Ward identity vanishes due to the agreement of A^ and A 2 . 
Second, rewriting Aa and A a ^ in the forms: 

Ai r) = A° + A* = V^c + gj&yfc (Al) 

A Hr) = A 10 + A ll = + g^Jr)^ (A2) 

(compare Fig.^| and Fig.|7|) we obtain ^Jn~ c [A\ — A l ^ r \ = ^/rr c [A 1 ^ — A l 2 °] = \q\n c and the 
third Ward identity reduces to 

"Xjnfau) = ^\Xjj(Q>v) +™n th ). (A3) 

Furthermore, in the first Ward identity the terms y/n^cu and (— l) v /n^a q 2 /2m appearing 
in the bracket on the right hand side are cancelled by the contributions u)A° a and |qr|/mA^° 
respectively. Since the contributions of T^Jjp are independent of a and (3 the difference 

— Egf* vanishes. Additionally we note that the differences of the Gross-Pitaevskii 
self-energies are given by Yf{ a — E° a = a\i resulting in a simplified expression of the first 
Ward identity 
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u;Ai 



q,uj) 



\q[ 

m 



(A4) 



Recalling the decomposition in Eqs. (|A"T, K2) and dividing A* and A^, 1 by their common 
factor gy/n~ c the proof of the first and second Ward identities is reduced to the check of the 
relation 



^Ann 

Using the decompositions (Fig. [| and |[ 



(r) = \Q[l(r) 



m 



Xnn ^nn 9Xnn) 

l(r) _ v i0 /■, , (r)\ 
Xjn Xjn ' 9Xnn) 



(A5) 



(A6) 
(A7) 



we just have to demonstrate the equivalent relation co>Xnn — (l < zl/ m ) Xj° n = 0> which can be 
done in a similar way as in |T2| (see also |TB[) by replacing the free-particle Green's functions 
used there by the Hartree-Fock Green's functions G HF of eq. Q2.16Q . In completing this proof 
we only need an energy dispersion law of the form e{p) ~ [p 2 /2m + const.} valid for the 
free-particle Green's function and the Hartree-Fock Green's function used by us. 



The third Ward identity can be simplified by the decompositions (Fig. ||] and ^) 



l(r) 
XjJ 

l(r) 
XjJ 



Xjn Xjn 



XjJ + Xjn' 



9Xmi. 



-Y 10 
Ann 
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a - gx 

which reduces its proof to the check of the condition 



X 1 °j 

\ A nJ 



(A8) 
(A9) 



ux l j°n(q,u) - ^x 1 jj{q^) = \q\ n t h- 
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m 



After multiplication with \q\/m we obtain: 



(A10) 



^xl(q^) ~ ^X l j°j(q,u; 



m 



m c 

q ( ,q 



fo(e(p)) - f (e{p + q)) 

m v 2 ' 



dp 3 
(27)3 

7^ii-(P + f ) I/o(e(p)) " /o(e(p + 9))] = -nth, 



— (P 

m 



2' 



which completes the proof. 
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FIG. 5. The chemical potential /U of the condensed phase in units of k^T as a function of 
the condensate fraction n c /n for fixed temperature. The plot is made for the coupling strength 
g = 0.3A 3 h £;£T/C(3/2). The horizontal long-dashed lines represent the borders of multi-stability, 
the upper one given by ^( max ) = 2 5 A t " h 3 C(3/2) the lower one // mm ) by the instability condition 
<?x(0, 0) = —1/2. The stable part of the curve is plotted as a solid line changing to a long-dashed 
line in the region of instability. 
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FIG. 6. The chemical potential \i of the condensed and the uncondensed Bose gas in units of 
ksT as a function of the particle density n in units of C(3/2)A^ 3 for the same value of the coupling 
constant g as in Fig. [|. The long-dashed curve represents the unstable part of the branch of the 
condensed phase. 
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FIG. 7. Diagrams contributing to A, 
The symbol ■ denotes the longitudinal component of the gradient. 
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FIG. 8. Regular current-density autocorrelation function Xj n corresponding to Xnn 
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FIG. 9. Regular current-current autocorrelation function Xjj corresponding to Xnn 
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FIG. 10. Real part (left) and imaginary part (right) of the different branches of the velocity 
of sound c = oj/\k\ from ( |3.1| ) and ( fO| ) in units y ksTc/m for the interaction gntot/ksTc = 0.1 
depending on the temperature T/T c . 
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FIG. 11. Real part and imaginary part of the different branches of the velocity of sound c in 
units -\JkBT c /m for the interaction gn tot /kBT c = 0.3 . 
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